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Abstract 



For the calculation of the partition function Z of small, isolated and 
interacting many body systems an improvement with respect to previous 
formulations is presented. By including inharmonic terms and employing 
CN ' a variational approach quantum effects can be treated even at very low 

temperatures. In order to test its accuracy the novel approach is applied 
, to the exactly solvable Lipkin-Meshkov-Glick model (LMGM). For ther- 

I modynamic properties and level densities good agreement with the exact 

^ ■ calculations is found. 

o 

■ 

1 Introduction 

Ch ' Level densities and thermodynamic properties of finite, interacting many body 

systems are of greater interest in many areas of physical research. A particular 
I challenge has been the temperature range where strong quantum effects render 

^ I mean field theories insufficient. Methods have been developed to approximate 

the fundamental partition function of systems with a separable two body inter- 
action. They are especially valuable as an arbitrary Hamiltonian with two body 
interaction V^'^^ can always be expanded into a series of separable terms of the 
form 

1 " 

n = f + = h^_j2k f^f" , (1) 

with one body operators T, H and F'^ (see e.g. |T|). Whereas the former must be 
hermitian, the effective two-body interaction might also be of the form {F^YF". 
In nuclear physics the quadrupole-quadrupole or the pairing interaction are ex- 
amples of this type. As one knows pair correlations also play a role for electronic 
systems, even if the latter are of mesoscopic scale as e.g. for ultra-small super 
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conducting metallic grains on which substantial experimental [2] and theoretical 
[31 E] research has focused recently. 

A very elegant approach to the calculation of the partition function is provided 
by the path integral formalism. Applying a Hubbard- Stratonovich transformation 
(see e.g. jH]) the two-body interaction may be turned into terms containing only 
one-body operators and scalar auxiliary fields. In the simplest approximation 
the latter may be treated at the mean field level, but extensions for the inclusion 
of many body quantum fluctuations are possible. They can be accounted for 
systematically in the sense of expansions about the semiclassical approximation. 
A static approximation was originally developed in order to study finite size 
effects in small superconductors 0. Later it had been applied in nuclear physics 
calculating thermodynamic properties and level densities of hot nuclei [HIH]. This 
method, which is valid in the high temperature limit, was coined the "Static Path 
Approximation" (SPA). Dynamical effects have then been included on the level 
of a local RPA making use of Gaussian approximations on top of the static 
path. Several version of this procedure exist under the names SPA -|- RPA [Hj, 
PSPA (for "Perturbed Static Path Approximation", which we take over in the 
following) P and Correlated Static Path Approximation (CSPA) [TniE]- These 
approaches are applicable in a temperature range where quantum effects are no 
longer negligible. 

However, at a certain temperature Tq the classical mean field solution may 
become unstable in functional space and large fluctuations render the Gaussian 
approximation insufficient. By thoroughly analyzing the appearance of instabili- 
ties for the many body problem at lower temperatures, it has been possible to go 
beyond the PSPA achieving in this way a smooth behavior around Tq. In this way 
the breakdown is shifted to the lower temperature To/2 In the present paper 
we want to report a novel approach [1^] which exploits a variational procedure. 
Its basic idea is taken over from ^1 where the case of a particle moving in 
a one dimensional potential has been treated. For the many body systems to 
be studied here, the method consists of two steps: First one goes beyond the 
Gaussian approximation by taking into account fourth order inharmonic terms. 
Second the local RPA modes are replaced by variational parameters. Their ap- 
propriate values are found by minimizing the free energy. Formally this method 
is applicable down to T = and delivers very good results for the free and inter- 
nal energy of the many body system. Unfortunately, at least for approximations 
adopted for our model the specific heat exhibits some unphysical features at very 
low temperatures; this problem will be addressed below in more detail. 

The paper is organized as follows: After reviewing the conventional approach 
to the partition function of interacting many body systems in Sect. |21we develop 
the novel method in detail in Sect. El In Sect. (Hit is tested at the exactly solvable 
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Lipkin-Meshkov-Glick model (LMGM) Finally the results are discussed and 
an outlook of future work is given. 

2 Partition function of an interacting many body 
system 

In the sequel we concentrate on the schematic special case of just one interaction 
term in and write 

n = H + ^FF. (2) 

Extension of the methods proposed in this paper to systems with more (indepen- 
dent) collective degrees of freedom {n > 1) is feasible in principle, although it 
may be tedious in detail. The coupling constant k = —\k\ is taken to be negative. 
In the nuclear case attractive interactions of this type generate collective modes 
of iso-scalar nature fTl]. In the present paper we will restrict ourselves to such 
a situation. With the path integral formulation repulsive interactions {k > 0) 
leading to iso-vector modes have been studied in [THI within the SPA and in [TH] 
on the basis of SPA+RPA. An extension of the methods proposed here to positive 
k should be possible extending the ideas of these references. A detailed study 
will be undertaken in a future publication. 

2.1 General formalism 

The partition function of the grand canonical ensemble reads 

Z{P) = Tr exp (-/3(7i - /.i)) , (3) 

where (3 = 1/T is the inverse temperature (in units with k^ = 1), n is the 
chemical potential and A the number operator. This 2{(3) shall be evaluated 
by functional integrals in imaginary time [Hj. As mentioned before, the form of 
the Hamiltonian in Q invites one to introduce the mean field approximation 
through a Hubbard-Stratonovich transformation (HST). Thereby the separable 
two body interaction FF is removed at the expense of introducing a collective 
variable through an auxiliary path ^(t). Since this procedure is well-known (see 
e.g. 011]) we simply state here the basic results which will serve as the starting 
point for our analysis. (We will stick to the notation used already in [201 with 
the exception of identifying C ^ ( .) 

After introducing the Fourier expansion g(r) = qo + J2r^o1r exp(iz/rr) with 
amplitudes obeying q_r = q* and the Matsubara frequencies Ur = {2n /hp)r = 
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(27iT/h)r the partition function may be written in a form containing a static 
part that depends on go only and a factor ({P, go) that corrects for effects of the 
dynamics: 



2n\k 



/ + 00 
dqo exp(-/3^^ff(/5,go)) (4) 



Here J-'^^{P, go) is understood as the effective free energy of the constituents of 
the total system (jSj) moving in a time dependent mean field g(r) of average go- 
The 

•^''^(Ago) = ^go'-^ln;^(/3,go) (5) 

is obtained in the limit where all time dependent fluctuations are neglected. 
Both quantities, go) and J-'^^^{P,qo), should not be mixed up with the 

free energy 

^(/5) = -llnZ(/5) (6) 

of the total self-bound system with Hamiltonian Q. In (jSI) there appears the 
grand canonical partition function z{P, go) which belongs to the static part 

nHST[qo]=H + qoF (7) 

of the mean field approximation to the Hamiltonian which comes in via the 
HST. Obviously, HkstIqo] is a sum of only one body operators; in the following 
its eigenenergies will be denoted by e;(go). 

All contributions from the dynamical part of the auxiliary path g(r) are con- 
tained in the factor C{f3,qo) which can formally be written as the functional 
integral 

C(Ago) = j V'q exp(-SE(/?,go)A) (8) 

(see e.g. P and [12]), where a slightly different notation is used. Please note that 
the Fourier representation of the measure T>'q basically consists of integrals over 
all qr but the go itself. To proceed further let us expand the Euclidean action se 
around the static path go: 




■^E = TTT I y]Ksqrqs+ 5^ prstqrqsqt+ (^rstuqrqsqtquj + 0{q^) (9) 

r,s,t^O r,s,t,ujtO / 
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Here the coefficients are given by ^3 



KsW,qo) = o'^'--^ + ^r^/ dTrdTs Lrs{(3,qo) , (10) 
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Prst{P,qo) = -^nr drrdrsdn Rrst{P,qo) , (H) 



o-rstu{(3,qo) = -jTn I dTrdTsdndTu [Srstuil3,qo) - SLrsiP,qo)LtuiP,qo)^l'^) 
with 

Lrs{f3,qo) = e'^-^e^'^-^ ^ (tF(r,)F(r,)),„ , (13) 

Rrsti/3,qo) = e'^-^e^'^-^e-'- 1 (tF(r,)F(r,)F(r,)),„ (14) 



and SrsUP,qo) = e'^^^^'^e^'^^^^e^'^^^'e^^"-" - (rF(r,,)F(r,)F(Ti)F(rJ),„ (15) 



1 

¥ 

Similar expressions have to be evaluated for the coefficients of higher order. In 
()13|1 to F{t) denotes the operator F in the interaction picture. T is the or- 
dering operator in imaginary time and (. . .)qg denotes the expectation value with 
respect to the canonical density operator for the Hamiltonian ((Tj). Notice that in 
this expansion no terms survive which involve the r-independent factor (F(r))q(,. 
Indeed, because of the presence of the exponential factors the r-integrations make 
such terms vanish. For this reason in Q terms linear in do not show up. The 
expectation values for products of several F{Tr) can be performed with the help 
of the Wick theorem. In this way these averages are traced back to contractions. 
After carrying out the r-integrations the expansion coefficients A^s, Prst and (Trstu 
can be expressed by one body Green's functions which are defined for the mean 
field Hamiltonian ((Zj). Details of this procedure can be found in [T^ IllSj. 



2.2 The conventional approach to quantum fluctuations 

As mentioned earlier, within the SPA all dynamical contributions are discarded, 
which formally is achieved by putting in (jH) the C^^^{P, qo) equal to unity jHllZlEl- 
For lower temperatures, when quantum properties tend to become important, 
fluctuations can be incorporated within the conventional version of the Perturbed 
Static Path Approximation (PSPA) jHl^^l], HI l^- There, the expansion 
is truncated after the second order terms in the g^, which effectively means 
to describe quantum effects on the level of local RPA. Using (|Tn|l and (fT!^ it 
can be shown that the following connection with the response function S{F)^ = 
—x{uj)Sq{uj) exists fIU\ : 

1 111 

Ks{P,qo) = 7T (l + ^X(i'^r))5r-s = -, 7^Sr-s = 7^ ^r{P , qo) ^r-s (16) 

2 2 1 - ftXcoll(l'^r) 2 
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In the second identity the general relation (see e.g. [T^) Xcou = + 

of the intrinsic to the collective response function has been used. Within the 

independent particle model (IPM) can be written as 

Uk>li^r + ^kM) 

Here, the ^^{(3, go) are the M local RPA frequencies which can be obtained from 
the secular equation [13 122] 

l + A;xK) = 0. (18) 

The huJki{qo) = efc(go) ~ are the frequencies of the intrinsic excitations and 
the primed product means that factors with = e; are omitted. Note that the 
numerator as well as the denominator of the product p7p consists of M factors. 
The largest contribution comes from the truly collective modes, where the 
differ essentially from the corresponding Uki- In PSPA all integrals to be evaluated 
in (jH)) are of Gaussian type, and cause no problem as long as Xr{P,qo) > for 
all r 7^ 0. For systems with instabilities one RPA frequency may become purely 
imaginary zuf < 0. Then the need for the Xr{P,qo) to be positive implies a 
condition for the temperature below which the conventional version of the PSPA 
breaks down as soon as the first one, Ai(/3, go), vanishes. This implies diverging gi- 
integrals due to large fluctuations in the gi-direction. In the theories of dissipative 
tunneling 1211 123 (for many body systems see [2011121113) the temperature 
To at which this breakdown happens is known as the "crossover temperature". 
For all T > To the PSPA is applicable with the dynamical factor C^^^^iPi qo) iii 
(0} being given by 

In fT2| the divergence of the gi-integrals at To has been cured by taking into 
account in the Se of ^ inharmonic terms up to fourth order in gi and q2- The 
method developed there and called "extended PSPA" (ePSPA) is very similar to 
the treatment of the crossover region in dissipative tunneling [23 • More precisely, 
in the expansion Q besides Xrs the three anharmonic coefficients pu-2, P-1-12 
and (Jii_i_i have been taken into account. The corresponding g2- and gi-integrals 
can be evaluated analytically, resulting in a lowering of the breakdown temper- 
ature down to To/2. At even lower temperatures the fluctuations in directions 
besides the gi become large, too. Due to the mutual coupling, in this case a 
similar, analytical treatment is no longer possible. In a further extension over 
the ePSPA, in ^21 a Low Temperature Approximation (LTA) has been proposed, 
that is applicable at any temperature. To reach this result the barrier region of 
J^^^^{P,qo) has been excluded from the go-integral in ([H). One deflciency of the 
LTA is that most observables show an unphysical sharp bend at T = Tq, which 
has its origin in this exclusion. 
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3 A variational approach 

Both the ePSPA and the LTA still have serious deficiencies in the low tempera- 
ture or the crossover region, respectively. Unfortunately, in the same temperature 
range physically interesting effects take place, like e.g. some known phase tran- 
sitions. For this reason a method is desirable which delivers reliable results also 
there. In P^IT^ a variational approach has been developed for the simpler system 
of a particle moving in a one dimensional potential, which allows one to evaluate 
the quantum partition function at arbitrary temperature to high accuracy. Using 
the formalism of coherent states, the underlying ideas have been applied to many 
body systems in Here, we follow a different approach and make direct use 
of the expansion Q. It is useful to rewrite the functional integral (jHI) for the 
dynamical correction by introducing a reference action sf^: 

C(/3,go) = Jv'q exp[-sf^/h] exp[-{sE-sf^)/h] 

= C(exp[-(.E-4'')M])^° (20) 

Later this reference action will be chosen such that the functional integral for the 
corresponding normalization factor = J Vq exp(— Sq/^) can be evaluated 
exactly. The average {■ ■ ■)q used in fEIHl is taken with respect to . Finally the 
inequality 

(exp[-(.E - sf;)/h])f; > exp[-(sE - sf^r^/h] (21) 

will be employed, which derives from the convexity of the exponential and can be 
taken over to functional integrals . The average on the right hand side of (PT|) 
is easier to calculate than the one on the left hand side. Therefore the strategy 
is to adjust a number of variational parameters in sf^ such that the right hand 
side is maximized. Using the correction factor 

C™(/3, go) = C eM-{sE - s^Y^/n] (22) 

this idea allows us to find an upper bound jF™''(/5, gg) > J^'^^(/3, qo) for the effective 
free energy defined in (jD) by minimizing 

r^'iP, go) = -^''"^(/5, go) - In Qo) (23) 

with respect to the variational parameters. The inequality 1)211) implies the fol- 
lowing restrictions to the partition function (j3)) and the free energy (jHI) of the 
total system: 



Z™^(/3) < Z{p) and J^™'"(/5) > J^ip) (24) 
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A convenient form for can be deduced from the PSPA action by the re- 
placement 



,PSPA 



hp 



— ^ Ar(/3,go) \qr 



r>0 



\k 



(25) 



r>0 



Here for the stiffnesses in g^-direction several choices are convenient. The two 
extreme cases are the following ones: 



A^(/3,go;^p 



k>l 



.K + UJ 



kV 



iyl + ^X) n 



(26) 



(b) 



In ) all M local RPA frequencies cjj^ appearing in (fTTj) are replaced by 

m = M adjustable parameters In case (b) only one adjustable parameter 



^If is introduced 



m 



1) which replaces that RPA mode zuf, which for systems 



with instabilities leads to the breakdown of the PSPA below Tq. This greatly 
simplifies the variational procedure by restricting the space spanned by the fi^ to 
a one dimensional one. It might also be possible to generalize to a small number 
1 < m <C M of RPA frequencies, namely those which belong to the modes with 
largest collectivity. 



For the normalization factor defined in fl2()|) the replacement 
leads to the following expression (comp. (fnijl ): 



r>0 



Ar 



(27) 



As we will see below, this form removes the convergence problems of the PSPA 
even in regions where the JF^^^ exhibits barriers. The correction factor (j^^ is 



evaluated best by decomposing the action se into se = +5se and calculating 



,PSPA 



Sq)q and (5se)q separately. With the definition 



the second order contribution is (see Appendix lA.lj) 



(2^ 



(b) 



^PSPA 



(29) 



r>0 
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The denominator of 11,, is two orders higher in z/^ ~ r than the numerator. 
Therefore for large r the terms fall off sufficiently fast for the sum to converge. 

In principle, when evaluating the {6se)q one needs to know the expansion 
coefficients in to all orders. In the case of a particle in a one dimensional 
potential this causes no problem as the coefficients can easily be obtained from 
the potential V{x) by mere differentiation For the many body problem on 

the other hand, more complicated expressions are needed, like e.g. (fTT|) and (fT^ 
with (fT^ -(fT3 j) . Due to symmetry arguments all contributions to {Sse)q vanish, 
which have an odd number of qr- To get a first orientation of how the novel 
variational approach works, we start from a truncated form s^'* = s^^^"^ + ^s^^ 
of the action se where only terms up to fourth order are taken into account. In 
Appendix IA.2l it is shown that the fourth order contribution to the average reads 

r,s>0 ^ * 

The behavior of the coefficients ars-r-s for large r and s is essential for the 
convergence of the sum (jHUj) . Details of the calculation of ^rstu are given in 
Sect. 4.1 of and in jj^]- It turns out that ars-r-s consists of a number of 
terms like 

111 

^rs-r-s ~ : : : ; (31) 

which fall off sufficiently fast to insure the convergence of (jHUI) . 

Putting together all contributions up to fourth order the correction factor (j^^ 
can be written in the form 

lnC(^)(Ago) = lnC-^(.^^^^-4°)??-^('54'^)?? (32) 

r>0 r>0 r,s>0 

This truncated form is used as an approximation to In (™ when minimizing jF™"" 
of (j2Sl)- The equations for the fi^ do not contain much physical information, 
for which reason they are not shown explicitly. It may be noted, however, that 
without any inharmonic terms {6se = 0) these equations are solved by the RPA 
frequencies fi^ = w"^ in the temperature range T > Tq. For the choice (b) in ()26p 
and ()28|) it is easy to show that = w\ indeed determines the only minimum 
of the corresponding JF™''. As expected, under these restrictions the variational 
procedure reproduces the PSPA. 
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4 Application to the Lipkin-Meshkov-Glick model 



There are only few models for many body systems which allow for an exact eval- 
uation of quantities of interest. One of them is the Lipkin-Meshkov-Glick model 
(LMGM) [in]. Using quasi-spin operators Jx and and a negative coupling 
constant k = — \k\ the Hamiltonian of the variant of the LMGM used here reads 

n = 2eJ, + 2k Jl . (33) 

Identifying H = 2eJz and F = 2J^ it has the structure of the Hamiltonian (j21) with 
one separable two body interaction. The corresponding mean field Hamiltonian 
dZj) reads ?tIhst[9o] = 2eJz + 2go^x- Its eigenvalues are supposed to be g-fo\d 
degenerate. For the our calculations we took g = 10 and e = 5 MeV. 

This model has often been used to test the results of SPA type approximations 
(see e.g. [HI IH HH 112] )• By varying the parameters of the model one gets a 
free energy J-'^^^{qo) which shows the desired property of developing a barrier 
at certain critical values, which may or may not be washed out by temperature 
smoothing. This is best seen after introducing the effective couphng constant 

^=\Ml>0_ (34) 
e 

As shown in Fig. ^ at zero temperature the J-'^^^{qo) of (jlj) has a barrier at 
go = for all K > 1 and none in the opposite case. The influence of temperature 
on this functional form is demonstrated by the dotted lines in Fig. |21 calculated 
for a K = 1.313, which was also used in [HIE]- The barrier disappears whenever 
(3 = 1/T becomes smaller than a critical value /3crit = 1/^crit- Another property 
of the LMGM is that it has only one RPA mode w. Under such conditions the 
choices (a) and (b) of and coincide and the optimization procedure of 
case (a) is facilitated greatly. 

Also shown in Fig. |21 are the corrections to the SPA which one gains by the 
approximations to the of (0J) described before. For these calculations the 
number of Fourier coefficients has been limited to = 100, and hence also 
for the products and sums in correction factors like (fT^ or (jH^ . It has been 
checked that larger values of A^ do not change the results anymore. Considerable 
corrections to the classical SPA are seen at any temperature for both the PSPA, 
the ePSPA as well as the variational approach. Up to (3 = 0.9 MeV~^ ~ 0.8/9o 
differences among them are hardly visible. The inverse crossover temperature for 
the PSPA lies at /3o = 1-139 MeV"^ At this temperature the sign change of the 
coefficient Ai takes place at the barrier top at go = of JF^^^ and leads to a 
divergence of the gi-integrals. For j3 > j3o this pathologic region grows to larger 
|go|. Deviations between the ePSPA and the variational approach are clearly 
visible at /? = 2 MeV~^ ^ 1.8/3o. 
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Figure 1: j^spa at T = as a function of go for different effective coupling 
constants k. 

It is remarkable that the effective free energy of the variational approach (fully 
drawn line in Fig.|21) shows only little structure in the barrier region of JF^^^. This 
feature does not change qualitatively at even larger (3 and is in agreement with 
the experience with the Feynman-Kleinert variational approach (FKV) fSJ |2Z] if 
applied to a one-dimensional double-well potential: There at low temperatures in 
the approximation Wi{xo) to the effective classical potential W{xo) (which can 
be seen as special one-dimensional cases of the effective free energies J-"™ and 
J^'^^, respectively) no barrier exists anymore. In j2HI Monte-Carlo simulations of 
the effective classical potential W{xo) have been discussed for the same problem. 
It is reassuring that at small temperatures also such a W{xq) was found flat in 
the barrier region. 

Next we turn to the free energy © of the total system with the partition 
function calculated from (j3)). The quality of the different approximations are 
illustrated in Fig. IHl From the left panel it becomes evident that for an effective 
coupling constant of k, = 1.313 the classical SPA is acceptable only at high tem- 
peratures (small P). In the regime P < 0.9 MeV^^ ~ 0.8/?o the PSPA improves 
this classical result by accounting for local RPA modes. A serious deficiency of 
the PSPA is its breakdown at p = po = 1.139 MeV"^ Even in a region below 
Po its results are no longer reliable. On the other hand the extension ePSPA of 
[T^ behaves completely regular in the crossover region /3 ~ /^q. It can be used 
as a reasonable approximation up to /3 ~ 1.2 MeV~^ ~ Po but breaks down at 
still larger p. In contrast, the variational approach is applicable even at very low 
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Figure 2: JF°^^ and different approximations to the effective free energy J-' 
of (0} for K = 1.313 as a function of go at various temperatures. The inverse 
critical and crossover temperature are /?crit = 0.4 MeV~^ and Po = 1.139 MeV~^ 
respectively. J^^^ is an even function of go- The arrows point to the smallest go 
where (at a given temperature) the PSPA still is formally applicable. 
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Figure 3: Left: Free energy of the total system in the LMGM as a function of (3. 
The arrow points to the inverse crossover temperature /3o- Right: Comparison of 
the accuracy of the PSPA and the variational approach to the free energy of the 
LMGM. 



temperatures. There the relative error is of the order of 1% only. 

Unfortunately, as can be seen from Fig. El the free energy !F^^\(3) of the vari- 
ational approach is not a monotonously increasing function of (3 as it should be 
according to the first law of thermodynamics. However, this deficiency is much 
worse for the PSPA and the ePSPA in a region T > Tq of formal applicability 
just before the breakdown. To the best of our knowledge this problem has never 
been payed attention to in the literature on PSPA. Conversely, the free energy 
of our variational approach still is an increasing function in the crossover region, 
where that of the PSPA and the ePSPA already decrease. As we will explain 
below, there are indications that the pathological feature of which oc- 

curs at even lower temperatures is due to the truncation of the reference action 
and not a property of the variational approach itself. Also, it is seen that for 
P > 1.2 MeV~^ ~ the variational free energy becomes smaller than the exact 
result. Although this behavior seems to be in contradiction to the restriction, 
in fact there is no violation of ()24|1 . To understand this, one should notice that 
the inequality (j2T| itself is valid for the truncated action s^^ as well. Of course, 
this inequality does not necessarily hold in case that different actions are used on 
both sides, say like se on the left and s^^ on the right. Hence, different to 
no inequality can be written down between the approximation J^'^^\(5) ~ J-'^^^{(3) 
and the exact free energy J^{P). 

We have been able to reproduce a similar behavior of the free energy using 
the FKV of ^3 l2Zj to calculate the partition function of a particle of mass M 
moving in a one-dimensional double-well potential V{x) = —ax^ -|- (1 — a)x^ 
(where < a < 1) with anharmonicity of higher than fourth order. Taking into 
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account the inharmonic terms to all orders the free energy of the FKV indeed 
is a monotonously increasing function of (3 and the inequality ()24|) is fulfilled 
by the variational approach for all temperatures. A truncation of the Euclidean 
action after the fourth order in the leads to a free energy that is smaller than 
the exact result and becomes decreasing at low temperatures. This observation 
indicates that also the unphysical decrease of the free energy J-'^'^\P) in the left 
panel of Fig. El does not reflect a deficiency of the variational approach as such, 
rather it is an artefact of the truncation of the reference action. 

In the right panel of Fig. El the accuracy of the PSPA and the variational 
approach is tested at the example of the free energy (jHl) by showing the ratio of the 
approximate results to the exact one. This is done for different effective coupling 
constants k > 1. It is observed that within its range of applicability, which is 
to say for (3 < Pq, the accuracy of the PSPA increases with increasing k. This 
may be traced back to the fact that the barrier height of the effective free energy 
J^^^^ at T = grows with k, see again Fig. ^ However, the larger the barrier the 
smaller the probability for quantum tunneling and, hence, the less important are 
fluctuations of large scale. However, when these fluctuations become more and 
more localized their treatment as local RPA modes becomes better and better. 
Contrary to the PSPA, the accuracy of the variational approach becomes worse 
the larger the barrier height. As can be seen from the right panel of Fig. El this 
feature is more pronounced for larger (3. In the range < (3 < jSo, on the other 
hand, the results of the variational approach are better than those of the PSPA. 
One should recall, perhaps, that nature of the inharmonic terms of J^^^^ taken 
into account by the variational approach still is based on a local approximation. 
It does not base on a large scale mode like the one necessary to describe quantum 
tunneling correctly. These findings are in agreement with the results of the FKV 
for the case of the one-dimensional double-well potential: There, it has been 
observed that for high barriers the FKV is unable to describe the splitting of low 
lying levels due to tunneling correctly [T3] . 

An interesting question is why the variational approach does not suffer a real 
breakdown. For an answer it is useful to look at the quantities 



They largely define the stiffnesses of s^^^^ and in gi-direction (see f03|) and 
mind Ui = 27t /h(3 as well as the fact that the product over /i reduces to a single 
factor in the LMGM). After the manipulations described in Sect. El and El f^es 
contributes to the denominators of the correction factors (fTn|l and (jH^ . Fig. HI 
shows the of (jHSj) as a function of P and for go = (where the square of 
the RPA-frequency zu'^ becomes negative for /? > Pent and takes on the largest 




with 




(35) 
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Figure 4: Comparison of the "stiffnesses" ilgg of Sg^^"^ and s^" in gi-direction 
defined in fl35|l using the RPA-frequency and the variational parameter Q, 
respectively. The arrow points to the inverse crossover temperature. The effective 
coupling strength is k = 1.313. 



absolute value). For = zu"^ the fl^f^ becomes negative at the inverse crossover 
temperature Po and the gi-integral needed in (jH)) diverges at larger p. Conse- 
quently, the formula (fTT^ . which has been derived under the condition that all 
g^-integrals converge, is no longer valid under these circumstances. On the other 
hand using = garanties fi^g to be strictly positive. As now all g^-integrals 
appearing in (^01) are of Gaussian type, the correction factor (jH^ of the varia- 
tional approach is well defined for all (3. As a result of this feature the effective 
free energy JF™'' shows much less structure than jF^^^"^, as already mentioned in 
connection with Fig. El 

Further tests of the accuracy of the different approximations may be done for 
the internal energy S{(3) = —d\aZ{(3)/df3 [13j of the total system or for the 
specific heat C{(3) = 9^ lnZ(/5)/9/3^. We concentrate on the latter for k > 1 
here. In Fig. El the classical SPA is seen to deliver a qualitatively correct behavior 
of the specific heat at all temperatures. For /3 ^ Pcrit it under-estimates the exact 
result by about 10%. At these smaller P, however, the PSPA and the ePSPA 
supply much better approximations. But they cease to be reliable already far 
below their formal limits of applicability at Pq and 2Pq respectively. There the 
variational approach leads to much better results, at least for not too large values 
of p. For P > 1.7 MeV ~ 1.5/3o, on the other hand, it leads to negative values. Let 
us note that this problem is much less pronounced for n < 1 and closely related 
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Figure 5: The specific heat of the LMGM as a function of /? in various approxima- 
tions. The arrows point to the inverse critical temperature Pent where the barrier 
in JF^P^ starts to build up at go = and the inverse crossover temperature Pq. 



to the appearance of positive curvatures of the J-'{(3) shown in the left panel of 
Fig. IHl The exact origin of this unphysical behavior is still unclear and must be 
investigated further. Perhaps and similar to the (apparent) violation of the 
effect might be related to the truncation of the expansion of the Euclidean 
action {Sse)q at fourth order in However, studying analog truncations in 

one-dimensional calculations within the FK V of fSl I2Z| , so far we have not found 
a similar behavior of the specific heat even for systems with strong inharmonic 
terms. Therefore the origin of the negative specific heat in Fig. El will have to be 
investigated further. 

One prominent application of the SPA or the PSPA has been the calculation 
of the level density p of finite, interacting many body systems beyond the in- 
dependent particle model [H El 13 In Fig- El this level density is shown as 
function of the excitation energy E* = E — Eq. The calculations have been per- 
formed with the Darwin-Fowler method (i.e. a saddle point approximation to 
the inverse Laplace transform from Z{f3) to p{E*), see e.g. p^). For the case 
of non-positive heat capacities this becomes problematic as then the equation 
E = = —d\iaZ{P)/dp no longer has unique solutions. In our apphcation 
we have concentrated on the one with smallest (3, simply because the other one 
lies in that regime of f3 where S{(3) has a positive slope, implying C{(3) < 0. For 
excitations E* above about 15 MeV (not shown in the plots) all approximations 
are in good agreement with the Darwin-Fowler calculation which is based on the 
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Figure 6: Logaritlimic representation of the Darwin-Fowler approximation to 
the level density as function of the excitation energy E* — E — Eq in various 
approximations. Whereas for k = 0.8 no barrier is present in J^^^^, for k = 1.313 
it builds up at /? = /?crit- 
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exact partition function Z{(3). At smaller E* , however, the level density of the 
classical SPA is qualitatively wrong. Here, all other approximations lead to con- 
siderable improvements over the SPA result. When the J-'^^^ has no barrier, like 
for K = 0.8 (see Fig. Q), the PSPA and ePSPA can hardly be distinguished from 
each other and show a qualitatively correct behavior. On the other hand, for 
K = 1.313, which is to say when the J^^^^ exhibits a barrier for low temperature 
(see Fig. both the PSPA as well the ePSPA do not reproduce correctly the 
behavior of the calculation using the full Z{j3). Nevertheless, the ePSPA shows 
a certain improvement over the PSPA. For n = 0.8 the variational approach re- 
produces perfectly well the calculation done for the exact Z{(3). For the case of 
a barrier with k = 1.313 one is in the situation discussed above where within 
our present approximation the specific heat becomes negative above a certain 
value of f3, see Figl^ Although a clear improvement over the ePSPA is seen, 
the present version of the variational approach is not capable of reproducing the 
qualitative behavior right at small E*. Finally we like to note that for E* below 
about 0.3 — 0.5 MeV the Darwin- Fowler method itself becomes questionable. This 
feature is in accord with the fact that the heat capacity approaches zero at about 
/? ~ 2 MeV-\ 



5 Discussion 



A variational approach to the partition function of a finite, interacting many 
body system has been developed and applied to the exactly solvable LMGM. It 
has been shown that for thermodynamic quantities like the free energy or the 
level density it delivers quite accurate results. This is especially important as 
the variational approach can formally be applied without any problems in the 
low temperature regime (3 > (3q, where the PSPA breaks down. In the crossover 
range /3 ~ /3o the variational approach considerably improves the approximations 
used so far. 

At very low temperatures there is still some deviation of the variational from 
the exact results, which may be explained as follows: Also the variational ap- 
proach starts from a static auxiliary field and takes into account small quantum 
fluctuations: g(r) ~ go- At very low temperatures, however, large scale fluctu- 
ations like quantum tunneling cannot be neglected. To include such effects the 
full nonlinearity of the problem must be considered. For interacting many body 
problems this is still an open problem. 

The range of applicability of the method proposed is not limited to problems 
of nuclear physics. Superconductivity in ultra-small metallic grains j2] is a prime 
example for possible applications in condensed matter physics, where up to now 
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the path integral approach has been hmited to the PSPA 

Unfortunately, the benefits gained through the variational approach must be 
payed for by a significantly larger computational effort. The two quantities that 
most influence the computation time are the number N of Fourier coefficients 
(where r = 1 . . . N) and the dimension m of the space spanned by the variational 
parameters (where fi = 1 . . .m and 1 < m < M) . With respect to the Fourier 
coefficients the computation time for the PSPA and the ePSPA of is limited by 
the calculation of the coefficients A,, of ()17|) (and the three coefficients Pu-2, 
P-i-12 of (fTT|l and aii_i_i of (fT2|l in the case of ePSPA). In contrast, for the 
variational approach in addition all the A^^ coefficients ars-r-s must be evaluated 
according to the formulas given in ^21 In the case of an extension of the 
expansion Q to order 2n even A^" coefficients would have to be calculated. Please 
have in mind that the one-dimensional analog jTHIT3ll7f] of the method proposed 
here does not suffer from this complication. There the expansion coefficients can 
be read off directly from the Taylor expansion of the underlying potential V{x). 

The computation time for the variational procedure required for minimizing 
with (p?^ is given by the task of finding the absolute minimum of a scalar 
function in an m-dimensional space. As this effort strongly depends on the under- 
lying function, no strict statement is possible concerning the number of necessary 
operations. Nevertheless it is intuitively clear that this number strongly increases 
with the dimension m of the corresponding space. Therefore, the choice (a) in 
()28|1 and may imply serious drawbacks. On the other hand, the choice (b) 
shows the way out of this problem. The essential idea is to restrict the variational 
procedure to the lowest-lying mode zuf that leads to the breakdown of the PSPA 
and contributes most to the product of (|T7j) . A similar philosophy is exploited 
in [29j in order to be able to introduce dissipative effects. In cases where the 
restriction of the number of the variational parameters to m = 1 is not good 
enough, it should be possible to introduce variational parameters fi^ for a small 
number m <^ M of such modes only and use the RPA frequencies cc^ for the rest. 
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A Evaluation of the averages 
A.l Second order contributions 



Using the definitions (j26|) - (j28|) the second order contribution to the average nec- 
essary on the right hand side of (j2Hl is 

/ PSPA go\go 



c 
hd 

^ J]n,A, (Re2(g,) + \^{qr)) 

r>0 
r>0 



|r / n ^ rfRe(g,)dIm(g,) exp ^- A ^ A, ■ {^^{q,) + Im2(g,)) 
1^1 

' ' r>0 

y ^ dRe(g,)dIm(g,) exp ("^^r ■ (Re2(g^) + ^(g,))^ (Re2(g,) + Im2(g,)) 
= /i^n^. (36) 

r>0 

Here in the second identity we have carried out all Gaussian integrals with s ^ 
r. In the third identity the remaining non-Gaussian integrals are evaluated. 
Note that the sum in the final expression is convergent due to the fact that the 
denominator of 11^ is two orders higher in Uj. = [2Ti/hj3)r than the numerator 
(see (P))). 
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A. 2 Fourth order contributions 



Due to symmetry arguments in fourth order the only contributing terms read 
= 7^ fUA\ dReiq,)dlm{q,) exp (-^ J] A, ■ (Re^(g,) + lm'{q,)) 



t>0 ' ' \ ' ' i>0 



•yft 
KB 

^ ars-r-s (Re^(gr) + Im^lgr)) {R^^ [qs] + Im^(g.)) 

T,S>0 



T,S>0 



' ' r,s>0 

y ':^\^^^^^^r)dl'ai{qr) exp ^-j^ A.^ ■ (Re^(gr) + Im^(g.r))^ (Re^(gr) + Im^(gr)) 
j -Ac?Re(g,)ciIm(g,) exp (-^A, ■ {Re\qs) + Im" {q,))^ (Re2(g,) + ^(g,)) 

- ^ 2^^— A, A/ ^"^^^ 

The quantities A,, defined in have been used. The special case r = s is 
included in (I37|l . Note that for constant coefficients cTrs-r-s the sum would not 
converge. Only the behavior (p?T|l garanties convergence of the sum. 
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